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We present a nonperturbative calculation of all multifractal scaling exponents at strong disorder 
for critical wavefunctions of Dirac fermions interacting with a non-Abelian random vector potential 
in two dimensions. The results, valid for an arbitrary number of fermionic flavours, are obtained by 
deriving from Conformal Field Theory an effective Gaussian model for the wavefunction amplitudes 
and mapping to the thermodynamics of a single particle in a random potential. Our spectrum 
confirms that the wavefunctions remain delocalized in the presence of strong disorder. 
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For a noninteracting disordered system near a metal- 
insulator transition, there now exists compelling evi- 
dence that wavefunctions exhibit multifractal behaviour 
at scales below the localization length L c . Namely, per- 
turbative renormalization-group analysis of replicated |l]] 
and sypersymmetric [|| nonlinear sigma models as well 
as numerical work (see for example |3[ ) , point to the exis- 
tence of irregular scaling properties of the local moments 
of the wavefunction distribution, 
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with respect to system size L < L c (a is the short- 
distance cutoff). The information regarding these irreg- 
ularities is most conveniently condensed into a scaling 
spectrum r(q) defined in terms of the above moments by 
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where the brackets denote disorder averaging over a suit- 
able ensemble, and Z(l) plays the role of the wavefunc- 
tion normalization. In this framework, a localized state 
has t(<7 > 0) = 0, while r(q) oc q for a simple fractal; any 
nonlinearity signals multifractality. 

There exists a well-studied class of exactly solvable toy 
models which are known to exhibit the phenomenon of 
multifractality, namely Dirac fermions subject to a ran- 
dom (Abelian or non-Abelian) vector potential (RVP) in 
2-d Euclidean space The Abelian case ||-||] cor- 

responds to a single flavour of 2-d Dirac fermions inter- 
acting with a random magnetic field, and applies to dis- 
ordered systems whose Fermi surface is effectively made 
of a single node characterized by a ^-shaped singular- 
ity at the Fermi energy (as for example the Chalker- 
Coddington network model, or degenerate semiconduc- 
tors Hj. Wavefunctions are known P,B[ to be localized 
in the vicinity of the Fermi energy except exactly at the 
Fermi level, where critical states exist. Their full multi- 
fractal spectrum was recently constructed in using 
a mapping to a generalized random energy model. 



On the other hand, the non-Abelian random vector po- 
tential (NARVP) model |p|-p^| describes systems with a 
Fermi surface collapsing to an arbitrary (even for a tight- 
binding Hamiltonian) number of Fermi points in the 2-d 
Brillouin zone. To each node then corresponds a differ- 
ent flavour of Dirac fermions. Besides being an effective 
theory for low-energy excitations of a d-wave supercon- 
ductor |[o|], the NARVP problem was recently obtained 
in the context of non-Hermitian Quantum Mechanics as 
the effective model for a 2-d particle in a random impu- 
rity potential subject to an imaginary (driving) vector 
potential H]. This last model has found applications 
in a wide variety of physical systems ranging from fluc- 
tuating vortex lines in superconductors with columnar 
defects [jl5[l to reaction diffusion problems in biological 
systems fig ] as well as advective diffusion in random me- 
dia JIt]]. Interestingly, numerical evidence [jl5| suggests 
the existence of a localization-delocalization transition of 
this 2-d system as a function of the driving potential, to 
be contrasted with Hermitian theories for which all wave- 
functions are expected to be localized in d < 2 |l8| . The 
nature of the critical states is however still an open issue. 

It is our purpose here to characterize the distribution of 
critical wavefunctions for the whole class of non-Abelian 
random vector potential models (defined on the group 
SU(N)), by achieving an exact calculation of their full 
multifractal scaling function r(q). We do this by reduc- 
ing the problem to the thermodynamics of a particle in a 
random potential used in the Abelian sector . At the 
fixed point of the theory located at the limit of infinitely 
strong disorder, we find that the multifractal spectrum 
r(q) for N flavours of Dirac fermions in a non-Abelian 
random vector potential is 
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for N = 2, 3, showing that zero-energy wavefunctions 
are always delocalized. This is (as was first alluded to 
in B) in contrast to the Abelian case |t]||, for which a 
sufficiently strong disorder localizes the wavefunctions. 
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In order to set the stage for the calculation in the non- 
Abelian case, let us briefly recall the procedure used in 
jj],D to solve the Abelian problem. The Hamiltonian 
for Dirac fermions interacting with an Abelian RVP is 
"Ha = ft + i A where the Pauli matrices are used as the 
2-d Euclidean Dirac matrices, i.e. A — = 1,2. 

The first step is to find an explicit expression for the 
zero-energy wavefunction for each realization of the disor- 
der. In the gauge where = e^d^cf), it can be checked 
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ensemble averaging on the RVP is made over a Gaussian 
white-noise distribution for the vector potential, i.e. 
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where g~A is the disorder strength. The key consists in 
interpreting Z(q) as a partition function for a single par- 
ticle in a random site potential V^(x) = 20(x), in other 
words to define the partition function and free energy 
as Z(/3) = X) x e" Wx) and F(/3) = -±lnZ(/3), where 
the sum is taken over the (L/a) 2 sites of the regularized 
system. By defining a microcanonical partition function 
(density of states) fl(E) counting the number of energy 
states in a window of width W around E 1 



n(E) = Y,5 w (E-V(x)), 
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where 5w(E) = e E l 2W , the partition function can 
be expressed in terms of the Laplace transform Z(J3) = 
J ^■Q,(E)e^' 3E . Since the disorder distribution (g) 
is a simple Gaussian, it is possible to calculate di- 
rectly the average (Cl(E)) = J T>(j>P[<f>]n(E) which de- 
pends solely on the two-point function (0(x)0(y)) = 

-^lnj |x ~^ 2 +a I where q c = ^2ir/g A . The num- 
ber of thermodynamic degrees of freedom being D — 
2ln(L/a), we can define the intensive width, energy and 
free energy as w = W/D,e = E/D,f (P) = F(j3)/D. 
The direct calculation M yields 
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We see that (O(e)) has the property of vanishing for 
|e| > e c = 2/q c in the thermodynamic limit D — > oo, 
and of diverging exponentially for |e| < e c . Moreover, 
it was explicitly shown in |^] that in the same limit, for 
a given realization of the disorder, the energy levels fall 
within the interval |e| < e c with probability 1. The en- 
tropy density was also argued to be selfaveraging, i.e. 
quenched and annealed averages coincide in this region. 

Since there are no states with e < — e c , the system 
undergoes a freezing transition. The free energy can be 
computed from basic thermodynamic relationships [p|, 



and the function r(q) is obtained using the correspon- 
dence r(q) = 2q\im a/L ^ [f Q (q) - / (1)]. In the weak 
disorder regime defined by q c > 1, it is given by 
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whereas in the strong disorder regime defined by q c < 1, 
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i.e. the parabolic law completely breaks down. This 
spectrum signals the localization of the wavefunction for 
strong enough disorder (r(q) vanishes for positive q). 

Let us now turn to the non- Abelian construction. The 
principles of our procedure will parallel the Abelian one, 
in needing to find explicit expressions for the zero-energy 
wavefunctions and an effective Gaussian-like field the- 
ory for the wavefunction modulus. This last step will 
be achieved by taking full advantage of the underlying 
group structure of the theory in order to extract only the 
elements essential to mult ifr act ality. 

We thus consider a Fermi surface whose low-energy ex- 
citations are effectively described by TV different flavours 
of Dirac fermions. The (27V x 27V) Hamiltonian is 
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where I is the N x N unit matrix acting in the fermion 
flavour space. The N x N matrices belong to the 
fundamental representation of the su(N) algebra (A^ = 
A^t a ,a — 1,..,7V 2 — 1), where t a are the generators of 
su(N) in the fundamental representation. We again use 
Pauli matrices as our 2-d Euclidean Dirac matrices, i.e. 
A = A^ (g) (T M ,/i = 1,2. It is convenient to separate the 
theory into chiral sectors, for which we define 

A±(x) = Ai(x)±ii4 2 (x) =ia ± « ?± (x) 5± 1 (x), (10) 

with the reality condition gZ 1 (x) = g+(x). The fields 
5±(x) are elements of the (noncompact) complexified 
group SU (TV) (obtained from SU (N) by letting the gen- 
erating coefficients take values in C instead of R). 

Zero energy eigenfunctions can be easily written down 
by decoupling the wavefunction from the vector fields 
with a non-Abelian gauge-like transformation: 
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where each of the two entries is composed of the ma- 
trix <?(x) multiplying TV-dimensional column vectors D 
whose components are entire functions of z, z respec- 
tively. These obey the zero eigenvalue Dirac equation 
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Note that this is just a convenient rewriting of the wave- 
functions introduced in [jl4j. Choosing constant {-D 1 * 1 }, 
the modulus of these eigenfunctions can be written as 
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where we have defined the new matrices h(x) = 
,g t (x) 3 (x). These belong to the coset SU C (N)/SU(N), 
as can be seen from their invariance under the left- 
multiplication of g(x) by u(x) € SU(N). 

The disorder averaging is here again made over a Gaus- 
sian white-noise distribution for the vector potential: 
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A striking fact about the NARVP model is that it is 
a theory which does not require the use of either the 
replica or SUSY approaches to perform the disorder aver- 
aging for zero-energy eigenfunctions |fl3|| . Moreover, the 
NARVP theory is known to have a conformally-invariant 
fixed point at infinite disorder strength |J. In this limit, 
the correct formulation of the disorder ensemble averag- 
ing leads to a functional integral over the coset space 
SU C (N)/SU(N) with a Wess-Zumino-Novikov-Witten 
(WZNW) action on level k = -2N for the field h(x) 
HfjJJ- The crucial point is that h(x), which is by con- 
struction a primary field of the coset WZNW model, is 
precisely the field that appears in jl3|). 

We will use the free field representation of this WZNW 
model, which can effectively be found from [^0|. Omit- 
ting the specifics, we here only state that the result is 
a theory for r free bosonic fields which we arrange in a 
vector 4> (r is the rank of the algebra, which is N — 1 in 
our case) together with residual fields. We find 



Sn = Tr - / d 2 x [N(d^, d^) + 2(p, <f,)H] + S Ie 
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where 1Z is the curvature (which can be interpreted as 
a background charge in the Dotsenko-Fateev construc- 
tion, and is specific to the non-Abelian problem), (.,.) 
represents the scalar product in the root space, and p is 
the Weyl vector (i.e. the half-sum of the positive roots 
of the algebra). S ICS is an action for the residual fields 
in the representation, which are uncoupled from the set 
4>. Since they will not influence the multifractality of the 
wavefunction, we can safely discard S ies . This free field 
representation has highest-weight primary fields p0|| 
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where {oij} is the set of simple roots, and A is the first 
fundamental weight, defined by the above properties. 

We know (see (|l]|)) that the modulus of the wave- 
function involves the matrix field /i(x), which is the pri- 
mary field of the WZNW model. We choose the sets 
{Z^} in such a way that the highest- weight field fffl ) 



appears in the modulus (this can always be done, since 
the WZNW model is invariant under rotation of the field 
by a constant group element). The problem is then 
greatly simplified by decomposing the vector </> in the 
(non-orthonormal) basis of the simple roots. The pri- 
mary field (wavefunction modulus) takes in that case the 
very convenient form of an exponential of only one com- 
ponent of the vector field 0, i.e. 
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The problem is that the different components of the 
vector field cj> are still dynamically coupled, since sim- 
ple roots are not orthogonal. Their scalar product is 



given by the Cartan matrix (ct^ay) = A 
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Since only the field <j)\ appears in the wavefunction mod- 
ulus, we wish to decouple it from the set {fa, i — 2, r}. 
This can be done by a simple nonsingular linear trans- 
formation (a Harish-Chandra decomposition), 



0(x) = P0'(x), 
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with Pij = 
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Thus, besides the unit diagonal, only the first column of 
P is taken to have nonvanishing elements. Moreover, <fii 
is left unchanged by the transformation (19). The matrix 
P T AP then becomes 
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with Cn — 2 — Pi 2 = jtzt- Using this transformation 
in (p^|), we achieve our main objective: the wavefunc- 
tion amplitude has completely decoupled from the set 
{</>,}, i = 2, r, which means that all information about 
multifractality is contained in the free field action (we 
have defined for convenience the new field (j> = —\<j>\) 



s N [4>] = 



d 2 x 



2n(N 



(21) 



representing the weight of the disorder distribution 
through P \<i>] tx e~ s W as far as the wavefunction modu- 
lus |* ( x )| 2 = e~ 2, * (x) /Z(l) is concerned. This is one of 
the central points of our paper. 

The action @ is a Gaussian Field Theory with curva- 
ture term, and can be viewed as the non-Abelian version 
of (|4|). We can thus perform a calculation for (Q(E)) 
similar to the one in M (see eq. (pp), for which we find 
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where c(x) = ^ J d 2 yy/g(y)TZ(y)G (y, x) and G = 
Go(x, x) is the cutoff limit of the regularized Green's 
function for the Laplacian. We have written c(x) as an in- 
tegral defined over a general manifold with metric (^(x), 
so the square root of the metric's determinant appears in 
the measure. For convenience, we take the topology of 
our system to be that of the sphere, and concentrate all 
the curvature on a single point that we send to infinity. In 
other words, we take yj g(x)72.(x) = 87rlimB^oo 5(B — x). 
The short-distance limit of the regularized Green's func- 
tion then becomes Go (x, y ) = — In j ^ x ~^ +a j . We 

find that the terms in ( p2| ) involving c(x) are subdomi- 
nant in the thermodynamic limit D — > oo, and that the 
proper expression for the disorder-averaged number of 
states in an energy interval ui around e (neglecting the 
width W as compared to terms of order ln(L/a)) becomes 
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with the critical energy now taking the value e c = 
y/2(N — 1)/N. This expression has exactly the same 
form as in the Abelian case || , and we can thus use pre- 
cisely the same thermodynamic-like reasoning to obtain 



r(q). We immediately see that q c = 2/e c 
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for the 



non-Abelian case, and that the full spectrum, in analogy 
with (|?]) , is given by ([3]) . q c is greater than one for every 
N, and we are thus always in the "weak disorder" regime 
of the Abelian solution, for which the logarithm of the 
wavefunction normalization factor Z(l) is selfaveraging 
and the wavefunctions are delocalized. Derealization of 
2-d Dirac fermions can also take place with other forms 
of disorder (see for a recent example [O] ) . 

The termination that we have obtained here is even 
more stringent than conjectured in fl2|j : there, arguments 
were given for a termination of r* (q) (calculated with 
it ^normalized wavefunctions) in general before q ~ N, 
while the case N = 2 was explicitly shown to have termi- 
nated between q — 2 and q = 3. Here, for N — 2, we find 
q c = 2V2, which is consistent with the results coming 
from the Logarithmic Conformal Field Theory Jl2| . 

In conclusion, we have studied the distribution char- 
acteristics of the critical wavefunctions of multiflavour 
Dirac fermions in a non-Abelian random vector poten- 
tial, by calculating exactly their full multifractal spec- 
trum for any number of flavours. This was achieved by 
building the relevant Gaussian field theory for wavefunc- 
tion amplitudes and using the mapping to the thermo- 
dynamics of a particle in a random potential as in the 
Abelian case. Our multifractal spectrum confirms that 
zero-energy wavefunctions are delocalized in the limit of 
infinite disorder strength, and yields further evidence for 



the existence of such types of delocalized states in 2-d 
non-Hermitian Quantum Mechanics. 
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